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A  physics-based  model  was  developed  for  nanosecond-pulse  discharges,  including  real¬ 
istic  air  kinetics,  electron  energy  transport,  and  compressible  bulk  gas  flow.  A  reduced 
plasma  kinetic  model  (23  species  and  50  processes)  was  developed  to  capture  the  dominant 
species  and  reactions  for  energy  storage  and  thermalization  in  the  discharge.  The  kinetic 
model  included  electronically  and  vibrationally  excited  species,  and  several  species  of  ions 
and  ground  state  neutrals.  The  governing  equations  included  the  Poisson  equation  for 
the  electric  potential,  diffusion  equations  for  each  neutral  species,  conservation  equations 
for  each  charged  species,  and  mass-averaged  conservation  equations  for  the  bulk  gas  flow. 

The  results  of  calculations  with  this  model  highlighted  the  path  of  energy  transfer  in  the 
discharge.  At  breakdown,  the  input  electrical  energy  was  transformed  over  a  timescale  on 
the  order  of  1  ns  into  chemical  energy  of  ions,  dissociation  products,  and  vibrationally  and 
electronically  excited  particles.  About  30%  of  this  energy  was  subsequently  thermalized 
over  a  timescale  of  10  /is.  Since  the  thermalization  time  scale  was  faster  than  the  acoustic 
time  scale,  the  heat  release  led  to  the  formation  of  weak  shock  waves  originating  near  the 
sheath  edge,  consistent  with  experimental  observations. 

I.  Introduction 

In  a  nanosecond-pulse  discharge,  the  input  energy  in  the  electrical  circuit  is  ultimately  converted  into 
heat  and  gas  motion.  Motivated  by  potential  uses  in  flow  control  applications,  we  are  interested  in  the  details 
of  the  physics  of  this  process. 

Repetitive,  pulsed  discharges  are  a  well-known  and  appealing  method  for  plasma  generation. 1-3  Such 
discharges  are  efficient  generators  of  both  ions  and  electronically  excited  species  because  of  their  high  instan¬ 
taneous  reduced  electric  field.4  In  aerospace  applications,  nanosecond-pulse  discharges  have  been  employed 
as  flow  control  actuators, 5-8  as  a  source  of  ionization  for  nonequilibrium  magnetohydrodynamic  devices,9,10 
and  as  a  means  for  enhancing  ignition  and  combustion.11, 12 

The  generation  of  shock  waves  by  volumetric  heat  release  in  pulsed  discharges  was  observed  and  explained 
in  the  1970s  in  the  context  of  gas  laser  technology.13, 14  Early  computations  by  Aleksandrov  et  al.14  assumed 
that  all  the  power  dissipated  in  the  discharge  immediately  went  into  heating  the  neutral  gas.  Popov,15 
however,  has  proposed  a  two-stage  heating  mechanism  in  which  product  species  and  electronically  excited 
species  are  generated  by  electron  impact,  and  then  the  stored  chemical  energy  is  converted  to  thermal  energy 
through  quenching  and  recombination  reactions.  (Reports  on  other  investigations  of  kinetic  mechanisms  can 
be  found  in  Refs.  16-18.)  Two-dimensional  calculations  have  been  carried  out  recently  by  Unfer  and  Boeuf,19 
assuming  instant  thermalization  of  30%  of  the  dissipated  power  that  goes  into  electronic  excitation. 

In  recent  experiments  by  Nishihara  et  al.,20  control  of  a  Mach  5  cylinder  flow  was  demonstrated  using 
a  pulsed  surface  dielectric  barrier  discharge.  The  hollow  cylinder  model  was  made  of  fused  quartz.  A  thin 
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copper  exposed  electrode  was  affixed  to  the  surface  of  the  cylinder,  with  a  second  copper  electrode  mounted 
inside.  A  combination  of  positive  and  negative  polarity  pulses  applied  to  the  two  electrodes  produced  a 
potential  difference  of  about  27  kV,  lasting  on  the  order  of  5  ns  (pulse  full- width  at  half  maximum).  The 
effects  of  the  energy  release  in  the  resulting  discharge  were  captured  using  phase-locked  schlieren  imaging. 
A  weak  shock  wave  was  seen  to  form  near  the  edge  of  the  exposed  electrode,  and  propagate  upstream  in  the 
shock-layer  flow  over  a  time  scale  on  the  order  of  microseconds.  When  the  perturbation  reached  the  bow 
shock,  its  shape  was  altered,  and  the  standoff  increased  by  up  to  25%. 

We  have  begun  to  formulate  a  high-fidelity  physical  model  of  the  energy  transfer  process  in  the  pulsed 
surface  dielectric  barrier  discharge.  For  simplicity,  we  have  focused  in  the  present  work  on  a  planar  geometry; 
experimental  evidence  shows  that  a  nearly  one-dimensional  discharge  can  occur  at  relatively  low  pressures.21 
Using  coupled  modeling  of  the  plasma  and  compressible  flow  in  a  one-dimensional  geometry  for  conditions 
representative  of  the  stagnation  region  of  the  Mach  5  cylinder  flow  experiments,  we  are  studying  the  dominant 
physical  effects,  including  energy  thermalization  kinetics  and  compression  wave  formation  and  propagation. 

To  this  end,  a  reduced  plasma  kinetic  model  (23  species  and  50  processes)  was  developed  first  by  carrying 
out  a  sensitivity  analysis  of  a  zero-dimensional  plasma  computation  with  an  extended  chemical  kinetic  model 
(46  species  and  395  processes).22  Transient,  one-dimensional  discharge  computations  were  then  carried  out 
using  the  reduced  kinetic  model,  incorporating  conservation  equations  for  each  species,  a  self-consistent 
computation  of  the  electric  potential  using  the  Poisson  equation,  and  a  mass-averaged  gas  dynamic  formu¬ 
lation  for  the  bulk  gas  motion.  This  paper  presents  the  major  results  of  these  calculations,  for  conditions 
corresponding  to  the  cylinder  flow  experiments. 

II.  Physical  Model 

This  section  presents  the  chemical  kinetics  model  and  the  conservation  laws  governing  each  species,  the 
bulk  gas,  and  the  electrodynamics.  There  are  two  main  differences  between  the  present  approach  and  other 
recent  studies  of  nanosecond  pulse  discharges  in  air.19, 23-26  First,  the  present  model  emphasizes  the  chemical 
kinetics  of  energy  storage  and  thermalization.  Second,  in  place  of  a  traditional  drift-diffusion  model  under 
the  local  field  approximation,  a  more  detailed  model  of  the  charged  particle  motion  was  employed  here  to 
try  to  capture  some  of  the  nonlocal  effects,  and  to  directly  address  the  transfer  of  momentum  and  energy 
between  the  charged  particles  and  neutrals. 

A.  Reduced  Kinetic  Model 

A  reduced  chemical  kinetics  model22  was  developed  to  identify  the  dominant  species  and  reactions  affecting 
the  energy  balance  and  the  rate  of  thermalization  in  the  discharge,  and  to  minimize  the  computational  cost 
of  the  transient,  one-dimensional  calculations  that  will  be  presented  here.  To  obtain  the  reduced  kinetic 
model,  we  applied  a  sensitivity  analysis  to  a  detailed,  transient,  zero-dimensional  air  plasma  model  used 
in  previous  work.2'  The  full  air  plasma  model  was  based  on  the  model  developed  by  Kossyi  et  al.28  The 
main  criterion  for  the  sensitivity  analysis  was  the  effect  of  individual  processes  on  the  time-dependent  energy 
fraction  thermalized  after  the  discharge  pulse.  The  species  included  in  the  model  are  given  in  Table  1  and 
the  reaction  mechanism  is  listed  in  Table  2. 

The  reduced  model  includes  the  following  species:  molecular  nitrogen  and  oxygen,  atomic  nitrogen  and 
oxygen,  ozone,  nitric  oxide,  electrons,  three  species  of  ions,  five  electronically  excited  neutral  species,  and 
eight  vibrationally  excited  neutrals.  The  reaction  mechanism  includes  electron  impact  reactions  (Reac¬ 
tions  1-10  and  43-50),  reactions  with  neutral  radicals  (Reactions  11-17),  quenching  reactions  and  other 
transitions  (Reactions  18-32),  charge  exchange  reactions  (Reactions  33-34),  and  electron- ion  recombination 
(Reactions  35-42).  Rate  coefficient  curve  fits  for  most  of  the  reactions  are  given  in  Table  2,  with  the  remaining 
coefficients  plotted  in  Fig.  1. 

For  the  high  reduced-electric-fields  considered  here,  the  rate  of  electron  impact  ionization  greatly  exceeds 
the  rate  of  electron  attachment,  and  processes  involving  negative  ions  do  not  affect  the  energy  balance.  Thus 
they  are  omitted  from  the  model. 

Further,  the  model  also  omits  vibrational  relaxation,  which  would  occur  over  time  scales  much  longer  than 
those  considered  in  the  simulations.  For  example,  for  a  temperature  of  300  K  and  a  pressure  of  0.05  atm,  the 
characteristic  time  scale  for  vibrational-translational  relaxation  of  vibrationally  excited  nitrogen  by  oxygen 
atoms29  is  about  200  /is.  The  time  scales  of  interest  in  the  present  simulations  are  more  than  an  order  of 
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magnitude  shorter  than  this. 

In  developing  the  model,  zero-dimensional  Boltzmann  equation  calculations  were  carried  out  with  a 
steady,  two-term  expansion  of  the  electron  energy  distribution  function  of  the  plasma  electrons.30  These 
calculations  employed  experimental  cross  sections  of  electron  impact  electronic  excitation,  dissociation,  ion¬ 
ization,  and  dissociative  attachment  processes.31,32  The  rate  coefficients  of  the  electron  impact  processes,  as 
functions  of  the  electron  temperature  Te,  were  derived  from  the  point  Boltzmann  solutions  by  averaging  the 
cross  sections  over  the  electron  energy  distribution  function.  The  results  of  these  calculations  are  summarized 
in  Fig.  1. 

Extremely  high  reduced  electric  fields  are  present  in  the  cathode  sheath.  Since  the  sheath  is  essentially 
particle- free,  the  critical  region  for  accurate  modeling  of  the  discharge  under  the  present  conditions  is  actually 
the  vicinity  of  the  cathode  sheath  edge,  where  significant  ionization  is  present  and  the  reduced  electric  field 
is  less  than  1000  Td. 

As  a  test  of  the  validity  of  the  model,  we  computed  the  Townsend  ionization  coefficient  a/N  =  ki/ue  for 
nitrogen  using  the  point  Boltzmann  solver,  and  compared  the  results  to  the  experimental  data  of  Haydon 
and  Williams.33  The  results  are  shown  in  Fig.  2.  The  solid  line  represents  the  results  of  the  Boltzmann 
calculations,  and  the  symbols  indicate  the  experimental  data.  Error  bars  are  included  on  the  experimental 
points,  representing  a  tolerance  of  three  times  the  experimental  standard  deviation.  Error  bars  are  omitted 
for  those  points  where  the  experimental  uncertainty  is  on  the  order  of  the  symbol  size. 

For  reduced  electric  fields  below  about  1000  Td,  the  calculations  and  measurements  agree  within  the 
experimental  uncertainty  (5%-17%).  For  fields  between  1000  Td  and  3000  Td,  the  Boltzmann  solution  is 
about  10%  higher  than  the  experimental  data.  Only  above  3000  Td  does  the  error  become  substantial. 
Thus  the  level  of  uncertainty  in  the  present  model  is  typical  of  continuum  electrical  discharge  computations. 
Moreover,  in  PIC-MCC  (particle-in-cell,  Monte-Carlo  collision)  computations  of  streamers  in  air,  Chanrion 
and  Neubert34  found  that  the  reduced  electric  field  had  to  exceed  1000  Td  for  runaway  electrons  to  be 
observed. 

Although  fields  higher  than  1000  Td  were  observed  in  the  present  calculations,  those  regions  were  essen¬ 
tially  free  of  charged  particles.  At  the  location  of  peak  charged  particle  number  density,  the  reduced  electric 
field  ranged  from  about  50  Td  at  5  kV  input  amplitude  to  130  Td  at  27  kV  input  amplitude.  Thus  the 
physical  model  for  ionization  at  the  sheath  edge  is  sufficiently  accurate  to  permit  basic  parametric  studies 
of  the  nanosecond-pulse,  dielectric  barrier  discharge. 


B.  Governing  Equations 

Rather  than  use  the  traditional  drift-diffusion  (one-moment)  model  for  charged  particle  motion,  here  we 
chose  to  employ  the  more  complex  five-moment  model35-40  in  order  capture  nonlocal  effects,  and  to  directly 
address  the  transfer  of  momentum  and  energy  between  the  charged  particles  and  neutrals. 

Mass,  momentum,  and  energy  conservation  equations  were  solved  for  each  charged  species: 

d 

—  (msns)  +  V-(msnsvs)  =  msws  (1) 

at 

d 

—  (msnsvs)  +  V-(msnsvsvs  +psI)  =  V-rs  +  qsns  E  +  As  (2) 

at 

d 

—  [msns£s\  +  V-[msnsvs£s  +psvs]  =  V-[rs-vs  -  Qs]  +  gsnsvs-E  +  Ms  (3) 

The  mass  per  particle  of  each  species  is  denoted  as  ms ,  and  the  corresponding  charge  per  particle  is  qs.  The 
species  number  density  is  ns,  the  velocity  is  vs,  the  total  energy  is  £s  =  es  +  ,  the  internal  energy  is  es, 

and  the  pressure  is  ps.  The  electric  field  is  E.  Here  ws  is  the  rate  of  production  of  particles  in  chemical 
reactions,  As  is  the  momentum  exchange  in  collisions,  and  Ms  is  the  energy  exchange.  The  species  shear 
stress  is  rs  and  the  species  heat  flux  is  Qs. 

The  pressure  is  found  from  ps  =  nskBTs ,  where  fee  is  the  Boltzmann  constant  and  Ts  is  the  translational 
temperature.  The  internal  energy  per  particle  is  assumed  to  have  the  form  mses  =  +  ks  (Ts  —  T0)/ (ys  —  1) , 
where  is  the  ratio  of  specific  heats  and  H0  is  the  heat  of  formation  per  particle  of  species-s.  Values  for 
the  heat  of  formation  of  each  species  are  given  in  Table  1. 

For  the  neutral  species,  a  diffusion  equation  formulation  was  used: 


dns 

~dt 


V- 


D, 


knTn 


SJps  =  uis 


(4) 
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where  Ds  is  the  diffusion  coefficient  and  w  is  the  mass-averaged  velocity. 

For  the  gas  as  a  whole,  a  mass-averaged  formulation  was  employed: 

v^+V-(pw)  =  0  (5) 

d 

7T :(pw)  +  V-(pww  +pl)  =  V-r  +  CE  (6) 

at 

Sr  (p£)  +  V-(pw£  +pw)  =  V- (r- w  —  Q)  +  Ej  (7) 

ot 

Here  p  is  the  overall  density,  p  is  the  overall  pressure,  (  is  the  space  charge,  and  j  is  the  electric  current.  The 
diffusion  velocity  for  each  species  is  Us  =  vs  —  w,  p£  =  pe  +  \pw2  +  J2S  |msns[/2  is  the  total  energy  of  the 
bulk  gas  flow,  and  e  is  the  total  internal  energy  per  unit  mass. 

The  formulation  allows  for  the  interaction  of  the  discharge  with  the  bulk  gas  flow  through  both  a  body 
force  and  heating.  Note  that  the  energy  delivered  to  the  bulk  gas  by  the  electric  field  is  initially  apportioned 
between  internal  energy,  thermal  energy,  and  kinetic  energy.  The  portion  that  goes  into  internal  energy  is 
set  by  the  chemical  kinetic  model,  which  also  allows  for  chemical  energy  to  thermalize  over  time. 

A  separate  temperature  was  computed  for  each  charged  species,  but  a  single  translational  temperature 
was  assigned  to  all  neutral  species.  The  neutral  translational  temperature  Tn  was  determined  by  subtracting 
the  flow  kinetic  energy  and  charged  particle  thermal  energy  from  the  total  energy  then  solving  the 
resulting  expression  for  neutral  gas  energy  for  the  translational  temperature. 

The  diffusive  fluxes  for  the  overall  gas  are: 


T  =  ^2  (T«  ~  ms«sUsUs)  (8) 

S 


The  summations  in  (8)-(9)  include  separate  flux  terms  for  each  charged  species,  but  the  fluxes  for  the  neutral 
species  were  combined  in  a  single  representative  term  for  air.  Note  that  we  have  retained  a  number  of  terms 
involving  the  diffusion  velocity  that  are  often  neglected  in  the  mass-averaged  formulation.  In  particular,  the 
electron  heat  conduction  and  energy  convected  by  the  electron  diffusion  flux  make  significant  contributions 
to  the  total  energy  balance  near  the  edge  of  the  cathode  sheath. 

The  Poisson  equation  was  employed  to  compute  the  electric  potential: 

V2<^  =  -C/co  (10) 


where  £q  is  the  permittivity  of  free  space.  The  electric  field  was  found  from  E  =  —  V</>. 


C.  Closure  Models 

In  the  equation  set  (1)— (3),  closure  models  are  needed  for  the  collision  source  terms,  the  viscous  stress  tensor, 
and  the  heat  flux  vector.  First  we  consider  the  elastic  and  inelastic  components  of  the  momentum  exchange 
term  As  =  Af  +  A(  and  the  energy  exchange  term  Ms  =  M f  +  M/. 

Because  the  gas  is  assumed  to  be  weakly  ionized,  the  primary  elastic  collisions  are  with  neutral  particles. 
The  following  models  are  used  for  the  elastic  components  of  the  collision  source  terms41,42  for  the  charged 
particles: 


Af  =  —nsmsvs(x s  -  v„)  (11) 

Mf  =  -ns  msVs  [3 kB{Ts  -  Tn)  +  (vs  -  vn)-(msvs  +  m„v„)]  (12) 

ms  +  mn 

Here  vs  is  the  rate  of  collision  with  the  neutrals,  and  mn  is  an  average  particle  mass  for  the  neutrals.  The 
neutral  translational  temperature  Tn  was  described  earlier,  and  vra  w  w  is  the  mass-averaged  velocity  of  the 
neutrals. 

Inelastic  momentum  and  energy  exchange  terms,  corresponding  to  the  reaction  mechanism  listed  in 
Table  2,  were  also  included  in  the  closure  model.  A  reaction-r  in  the  kinetic  model  can  be  expressed 
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symbolically  in  the  form: 


N  N 

a'rsMs  — >  ^  v"s Ma  (13) 

s=l  s=l 

where  v'rs  and  v"s  are  the  stoichiometric  coefficients  for  the  reactants  and  products  and  M.s  is  the  chemical 
symbol  for  species-s.  Using  the  law  of  mass  action,43  the  rate  of  production  of  species-s  can  be  expressed  as: 

N 

Us  =  iVrs  ~  Ks  )  kr  ^  ( 14) 

r  i=  1 


where  kr  is  the  reaction-rate  constant. 

The  form  for  the  inelastic  momentum  source  terms  was  A4  =  mswsvs  for  all  species,  and  the  energy 
source  term  for  the  ions  was  Ml  =  msus£s.  A  more  complex  energy  exchange  term  was  necessary  for  the 
electrons: 

Ml  =  msujs£s  —  RrAHrs  (15) 

r 

Here  Rr  =  kr  JlfLi  n«rs  is  the  rate  of  progress  of  reaction-r  and  A Hrs  is  the  corresponding  energy  loss  for 
species-s.  The  energy  loss  includes  the  cost  of  heating  the  electrons  to  the  average  electron  temperature, 
and  is  given  for  each  reaction  in  Table  3. 

For  a  plasma  in  a  uniform  electric  held,  as  considered  in  the  Boltzmann  solutions  used  to  establish  the 
rates  of  the  electron  impact  reactions,  overall  energy  conservation  reduces  to  balance  between  the  work  done 
by  the  electric  held  and  the  energy  exchange  in  collisions.  This  balance  was  checked  for  consistency  for  both 
the  electron  (3)  and  total  (7)  energy  equations,  and  an  error  of  no  more  than  a  few  percent  was  observed 
over  a  range  of  electron  temperatures  up  to  100  eV. 

For  the  hux  terms,  it  was  assumed  that  the  viscous  term  had  a  Newtonian  form,  with  Stokes  hypothesis 
applied,  and  that  the  heat  hux  followed  Fourier’s  law: 


1~s  —  (Ay 


(Vvs)  +  (Vvs)T-=V-vs  I 
Qs  =  - ksX7Ts 


where  /ivs  is  the  viscosity  and  ks  is  the  thermal  conductivity  for  species-s. 


(16) 

(17) 


D.  Transport  Properties 

The  viscosity  and  thermal  conductivity  of  the  neutrals  were  found  from  standard  correlations  for  air.44  The 
ion-neutral  collision  rate  was  determined  as  follows:45 


vs/N  =  crs 


8  kBTs 


(18) 


where  crs  is  an  effective  collision  cross-section  and  N  is  the  overall  number  density.  The  collision  cross-sections 
used  for  each  species  are  given  in  Table  1;  the  values  were  determined  from  mobility  and  diffusion  coefficient 
data  in  the  literature.45-51  As  discussed  previously,  the  collision  rate  for  the  electrons  was  determined  from 
a  point  Boltzmann  solution  as  a  function  of  electron  temperature  ve/N  =  f(Te),  and  is  shown  in  Fig.  1. 

The  relationship  between  the  collision  rate  and  the  transport  coefficients  is  as  follows: 


Ds 


kBTs 

ms^s 


l^vs  —  ^nsmsDs 

5 

ks  =  -kBnsDs 


(19) 

(20) 
(21) 


The  forms  for  the  viscosity  and  thermal  conductivity  correspond  to  a  Lewis  number  of  unity  and  a  Prandtl 
number  of  2/3.  Note  that  the  viscosity  and  thermal  conductivity  coefficients  are  proportional  to  number 
density.  These  transport  properties  are  small  for  the  low  charged  particle  number  densities  in  the  cathode 
sheath,  but  become  significant  near  the  sheath  edge. 
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III.  Numerical  Implementation 


This  section  describes  the  numerical  algorithm  and  boundary  conditions.  In  short,  an  implicit,  second- 
order  upwind  scheme  was  employed,  along  with  standard  boundary  conditions  for  the  dielectric  surfaces. 

A.  Numerical  Methods 

The  calculations  were  carried  out  using  the  Air  Force  Research  Laboratory  computer  code  HOPS  (Higher 
Order  Plasma  Solver),  which  includes  several  physical  models  and  numerical  schemes.22,52  Here,  the  physical 
model  consisting  of  (l)-(7)  and  (10)  was  solved  using  an  implicit,  second-order,  upwind  formulation.  The 
equations  were  solved  in  a  nondimensional  form  that  has  been  described  in  previous  papers.53 

Time  integration  of  the  conservation  equations  (l)-(7)  was  carried  out  using  a  second-order  implicit 
scheme,  based  on  a  three-point  backward  difference  of  the  time  terms.  The  formulation  is  similar  to  the 
standard  technique  of  Beam  and  Warming,54  but  is  adapted  here  to  a  multi-fluid  formulation  with  different 
models  for  particle  motion.  The  solution  was  time  accurate,  with  all  equations  advanced  with  the  same  time 
step. 

Approximate  factoring  and  quasi-Newton  subiterations  were  employed.  The  implicit  terms  were  linearized 
in  the  standard  ‘thin  layer’  manner.  The  implicit  terms  were  evaluated  with  second-order  spatial  accuracy, 
yielding  a  block  tridiagonal  system  of  equations  for  each  factor.  The  species  were  loosely  coupled,  limiting 
the  rank  of  the  flux  Jacobian  matrices  to  the  order  of  the  moment  model  (one  for  the  diffusion  equation 
formulation,  five  for  the  charged  particles  and  overall  conservation  equations).  The  resulting  equations  were 
solved  using  a  standard  block  tridiagonal  solver,  and  the  change  in  the  solution  vector  of  conserved  variables 
was  driven  to  zero  by  the  subiteration  procedure  at  each  time  step.  Three  applications  of  the  flow  solver  per 
time-step  were  employed  for  the  present  work. 

For  the  charged  particles  and  the  bulk  gas,  the  Roe  scheme55-57  was  employed  for  the  inviscid  fluxes.  For 
the  diffusion  model,  a  simple  upwinding  scheme  was  employed,  based  on  the  convection  velocity.  In  both 
formulations,  stability  was  enforced  using  the  minmod  limiter  in  the  MUSCL  formalism.58 

The  Poisson  equation  (10)  was  solved  at  the  end  of  each  subiteration  in  the  implicit  time-marching 
scheme.  (This  procedure  yields  a  stable  time  step  that  is  comparable  to  that  obtained  using  methods  based 
on  the  linearization  of  the  right-hand-side  of  the  Poisson  equation.59)  The  numerical  scheme  was  adapted 
from  the  approach  described  by  Holst.60  The  formulation  of  the  implicit  scheme  was  analogous  to  that 
used  for  the  conservation  equations,  with  linearization  of  the  implicit  terms,  approximate  factoring,  and  an 
iterative  procedure  that  drives  the  change  in  the  solution  to  zero.  The  spatial  derivatives  were  evaluated 
using  second-order  central  differences,  and  the  system  was  solved  using  a  standard  tridiagonal  algorithm. 

B.  Boundary  Conditions 

A  one-dimensional  computational  domain  was  employed  to  represent  a  double  dielectric  barrier  discharge 
(Fig.  3).  All  the  boundary  conditions  were  enforced  to  second-order  numerical  accuracy.  No-slip  boundary 
conditions  with  a  constant  temperature  wall  were  employed  for  the  bulk  gas.  A  zero  wall-normal  derivative 
was  imposed  for  the  neutral  species. 

Standard  boundary  conditions  were  employed  for  the  charged  particles.  First,  provisional  conditions  at 
the  wall  were  determined  by  setting  the  normal  derivative  to  zero.  Then,  if  the  ion  flow  was  away  from  the 
boundary,  the  ion  flux  was  set  to  zero  and  the  ion  temperature  was  set  to  the  wall  temperature.  For  electron 
flow  from  the  wall  into  the  domain,  the  normal  component  of  the  electron  velocity  was  set  to  that  required  to 
satisfy  secondary  emission,  and  the  electron  temperature  was  set  to  a  specified  emission  temperature  (here, 
1  eV). 

Simplified  boundary  conditions  were  employed  to  model  a  thin  dielectric  electrode  coating.49  The  dielec¬ 
tric  layer  was  assumed  to  be  sufficiently  thin  to  be  approximated  a  linear  potential  profile  (uniform  electric 
field  Ed).  The  electric  field  inside  the  dielectric  was  related  to  the  electric  field  E  at  the  surface  through  the 
relation  eoE  —  ereoEd  =  ern,  where  a  is  the  surface  charge  density  and  n  is  the  unit  normal  vector  pointing 
into  the  computational  domain.  The  surface  charge  was  determined  by  integrating  da / dt  =  —  j  •  n  for  each 
surface  point,61  using  a  time-marching  scheme  analogous  to  that  of  the  main  governing  equations. 
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C.  Two-Stage  Formulation 

For  computational  efficiency,  the  calculations  were  carried  out  in  two  stages.  The  first  stage  encompassed 
the  first  400  ns  of  the  discharge.  For  this  stage,  the  full  physical  model  discussed  above  was  employed. 
Since  electromagnetic  effects  and  charged  particle  motion  become  negligible  after  the  input  pulse  dies  away 
(after  about  24  ns),  in  the  second  stage  of  the  computations  (0.4  /is  to  100.4  /ns),  the  electric  field  was  set 
to  zero  and  neutrality  was  enforced  by  appropriately  setting  the  electron  number  density.  This  two-stage 
approach  resulted  in  a  substantial  savings  in  computational  cost,  and  tests  comparing  this  approach  to  the 
full  formulation  showed  a  negligible  difference. 

IV.  Results 

Calculations  were  carried  out  for  a  one-dimensional  discharge  under  conditions  representative  of  the  stag¬ 
nation  region  of  a  Mach  5  cylinder  flow  experiment  described  in  a  previous  publication.20  The  corresponding 
stagnation  conditions  were  4.74  kPa  (36  Torr)  and  310.3  K;  for  each  case  in  the  calculations  discussed  below, 
the  initial,  uniform  state  of  the  neutral  gas  was  set  to  these  values. 

The  configuration  considered  here  is  illustrated  in  Fig.  3.  The  problem  is  one-dimensional.  In  the 
simulations,  the  right  electrode  was  grounded,  and  the  left  electrode  was  powered  with  the  input  signal 
Vs  =  —Vo  exp  [— (f  —  to)2/r2] .  For  the  present  calculations,  the  input  signal  parameters  were  Vo  =  5  kV  to 
27  kV,  r  =  3  ns,  and  to  =  12  ns. 

Both  electrodes  were  assumed  to  be  covered  with  a  d  =  1  mm  thick  dielectric  coating,  with  a  relative 
dielectric  constant  of  er  =  3.8,  chosen  to  be  representative  of  fused  quartz.  The  discharge  gap  (distance 
between  the  dielectric  surfaces)  was  taken  to  be  L  =  2  mm.  The  secondary  emission  coefficient  was  7sem  = 
0.05. 

The  initial  mole  fraction  of  the  electrons  and  each  of  the  neutral  minor  species  was  taken  to  be  1  x  10-10. 
The  number  densities  of  the  vibrational  states  of  N2  were  set  to  an  equilibrium  distribution  for  the  initial 
temperature.  The  mole  fraction  for  each  ion  species  was  equal,  and  set  so  that  the  space  charge  was  zero. 
The  initial  electric  field  was  zero.  The  ion,  electron,  and  neutral  temperatures  were  set  to  310.3  K,  and  all 
velocities  were  set  to  zero. 

For  the  calculations  presented  here,  a  uniform  grid  of  1001  points  across  the  gap  was  employed  (Ax  = 
2  /im).  Grid  resolution  studies  presented  in  a  previous  paper22  indicate  that  this  level  of  grid  resolution  is 
sufficient  for  this  problem. 

The  time  step  for  the  Stage  1  calculations  (see  Sec.  III.C)  was  taken  to  be  0.5  ps  for  the  5  kV  case. 
The  minimum  chemical  kinetic  time  scale  and  dielectric  relaxation  time  scale  were  both  about  20  ps,  and 
the  minimum  electron  diffusion  time  scale  was  about  2.5  ps,  so  the  time-evolution  of  the  discharge  was 
well-resolved  in  the  calculations.  The  corresponding  time  step  for  the  10  kV  and  20  kV  cases  was  0.25  ps, 
and  for  the  27  kV  case  it  was  0.125  ps.  For  the  Stage  2  calculations,  the  time  step  was  taken  to  be  0.5  ns 
for  all  cases. 

A.  Baseline  Case 

A  baseline  case  with  Vo  =  5  kV  was  studied  first.  The  general  structure  of  the  discharge  is  shown  in  figures  4- 
5,  which  show  profiles  of  number  density,  potential,  velocity,  and  temperature  at  times  corresponding  to  the 
peak  in  the  input  voltage  (12  ns,  Fig.  4)  and  the  effective  end  of  the  input  signal  (24  ns,  Fig.  5).  The  cathode 
sheath  is  evident  at  the  left  in  the  figures  as  a  region  largely  free  of  charged  particles.  There  is  a  relatively 
large  electric  field  in  the  sheath,  resulting  from  the  region  of  space  charge  at  the  sheath  edge  (x/L  ss  0.3). 

Note  that  the  electron  velocity  and  temperature  at  the  left  boundary  are  set  by  secondary  emission.  At 
the  peak  of  the  input  signal  (12  ns),  ion  velocities  reach  ~  103  m/s  in  the  sheath,  and  electron  velocities 
~  106  m/s.  Corresponding  temperatures  reach  ~  1  eV  for  the  ions  and  50  eV  for  the  electrons.  Although 
these  values  are  quite  high,  the  charged  particle  density  is  extremely  low  in  the  left  sheath,  and  there  is 
essentially  no  conduction  current  or  energy  storage  there.  Despite  the  high  dimensional  velocities,  it  should 
be  noted  that  the  maximum  electron  Mach  number  is  about  0.7,  and  the  maximum  ion  Mach  number  lies 
in  the  range  of  1.0-1. 6. 

As  the  input  signal  dies  away,  the  discharge  enters  a  dual-cathode  regime,4  as  seen  in  Fig.  5  at  24  ns.  At 
both  boundaries,  ions  travel  toward  the  wall,  and  electrons  away.  At  this  stage,  the  space  charge  in  the  left 
sheath  has  decayed  significantly,  but  the  right  sheath  is  quite  prominent. 
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Note  the  qualitative  differences  between  the  left  and  right  sheaths.  The  left  sheath  is  formed  during  the 
input  voltage  rise  by  way  of  an  ionization  wave  propagating  into  a  very  weakly  pre-ionized  plasma.  The  right 
sheath  forms  during  the  voltage  fall,  in  a  region  with  significant  existing  ionization.  Its  characteristics  are 
similar  to  those  of  sheaths  at  a  negative  electrode  in  a  DC  or  RF  discharge.  In  particular,  there  is  significant 
ion  current  in  the  right  sheath,  in  contrast  to  the  absence  of  conduction  current  in  the  left  sheath. 

B.  Acoustic  Waves 

The  interaction  between  the  discharge  and  the  bulk  gas  flow  occurs  through  the  source  terms  in  Eqs.  (6)-(7). 
Terms  of  interest  include  the  body  force  £E,  the  rate  of  work  done  by  the  body  force  </E-w,  and  the  total 
electromagnetic  power  Ej.  Maximum  values  of  these  terms  occur  near  the  edge  of  the  left  sheath  at  the  peak 
of  the  input  waveform,  and  near  the  edge  of  the  right  sheath  as  the  input  waveform  decays.  The  body  force 
is  always  directed  toward  the  electrodes,  and  tends  to  be  stronger  for  the  left  sheath.  For  the  Vo  =  5  kV 
input  amplitude,  the  order  of  magnitude  of  the  terms  was  as  follows:  peak  body  force  105  N/m3,  rate  of 
work  done  by  body  force  104  W/m3,  and  total  electromagnetic  power  1010  W/m3.  The  corresponding  values 
for  Vo  =  27  kV  were:  peak  body  force  107  N/m3,  rate  of  work  done  by  body  force  10 7  W/m3,  and  total 
electromagnetic  power  1011  W/m3.  (These  large  values  exist  only  for  time  scales  on  the  order  of  the  width 
of  the  input  waveform,  and  for  spatial  scales  on  the  order  of  the  discharge  volume.)  Since  the  total  power  is 
four  to  six  orders  of  magnitude  larger  than  the  rate  of  work  done,  acoustic  waves  that  appear  in  the  discharge 
must  primarily  be  a  result  of  thermal  energy  release. 

The  formation  of  acoustic  waves  from  the  heat  release  is  illustrated  in  Fig.  6,  which  shows  temperature 
and  velocity  profiles  at  several  times  for  the  5  kV  and  27  kV  cases.  In  both  cases,  heating  at  the  cathode 
sheath  edge  (x/L  ss  0.3)  generates  waves  that  propagate  away  from  that  location.  Similarly,  heating  near 
the  right  boundary  creates  a  wave  that  travels  into  the  domain. 

Focusing  on  the  vicinity  of  the  left  sheath  edge  at  0.4  fis  for  the  5  kV  case  (red  curve  in  Fig.  6b),  we  note 
an  asymmetry  in  the  waveform  that  is  probably  a  result  of  the  left-directed  body  force.  The  velocity  of  the 
left-running  wave  is  about  60%  higher  than  the  that  of  the  right-running  wave.  The  corresponding  difference 
for  the  27  kV  case  (red  curve  in  Fig.  6d)  is  about  40%.  Nonetheless,  the  body  force  has  a  relatively  small 
influence  on  the  formation  of  the  acoustic  waves. 

The  results  are  qualitatively  similar  for  all  the  cases  examined  here,  with  the  wave  amplitude  and  speed 
increasing  with  the  amplitude  of  the  input  signal.  The  wave  speeds  are  nearly  sonic  for  all  cases;  the  Mach 
number  does  not  exceed  1.08.  (The  sound  speed  is  353  m/s  for  air  at  310.3  K,  and  the  wave  speeds  lie  in 
the  range  of  350-380  m/s.)  At  later  stages  in  the  computations  (not  shown),  the  waves  reflect  from  the  solid 
boundaries,  producing  complex  wave  patterns. 

These  computational  results  are  consistent  with  the  waves  emanating  from  the  cathode  and  anode  in  a 
pulsed  CO2  laser  observed  by  Pugh  et  al.13  using  interferometry.  Those  authors  noted  both  generation  of 
waves  by  discharge  heating  and  reflection  of  the  resulting  waves  from  solid  boundaries.  The  configuration  in 
the  Mach  5  cylinder  flow  experiments20  was  somewhat  more  complicated,  but  similar  results  were  obtained. 
A  weak  shock  formed  near  the  edge  of  the  exposed  electrode  and  propagated  into  the  shock  layer  flow. 
The  results  of  the  present  calculations  are  qualitatively  consistent  with  those  of  both  experimental  studies. 
Further,  wave  speeds  extracted  from  the  computational  results  shown  in  Fig.  6d  lie  in  the  range  of  350- 
380  m/s,  in  general  agreement  with  the  value  of  375  m/s  measured  in  the  Mach  5  cylinder  flow  experiments. 

V.  Summary  and  Conclusions 

Numerical  calculations  were  carried  out  to  examine  the  physics  of  the  operation  of  a  nanosecond-pulse, 
dielectric  barrier  discharge  in  a  configuration  with  planar  symmetry.  This  simplified  configuration  was  chosen 
as  a  vehicle  to  develop  a  physics-based  nanosecond  discharge  model,  including  realistic  air  plasma  chemistry 
and  compressible  bulk  gas  flow.  Discharge  parameters  (temperature,  pressure,  and  input  waveform)  were 
selected  to  be  representative  of  recent  experiments  on  bow  shock  control  with  a  nanosecond  discharge  in  a 
Mach  5  cylinder  flow. 

The  computational  results  qualitatively  reproduce  many  of  the  features  observed  in  the  experiments, 
including  the  rapid  thermalization  of  the  input  electrical  energy  and  the  consequent  formation  of  weak  shock 
waves.  Despite  the  simplifications  employed  in  the  calculations,  the  computed  translational  temperature 
rise  (40  K)  and  nitrogen  vibrational  temperature  rise  (370  K)  were  of  the  same  order  of  magnitude  as  those 
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measured  experimentally  (50  K  and  500  K,  respectively).  The  results  illustrate  how  input  electrical  energy 
is  rapidly  transformed  (over  roughly  1  ns)  at  breakdown  into  ionization  products,  dissociation  products,  and 
electronically  excited  particles,  and  how  thermalization  occurs  over  a  relatively  longer  time-scale  (roughly 
10  fi s).  The  quenching  of  electronically  excited  states  and  electron-ion  recombination  make  roughly  equal 
contributions  to  the  thermalization  of  stored  chemical  energy. 

This  work  represents  a  first  step  towards  detailed  modeling  of  the  nanosecond-pulse  dielectric  barrier 
discharge  actuator.  This  type  of  actuator  shows  promise  for  high-speed  flow  control,  and  accurate  numerical 
modeling  will  contribute  the  optimization  of  these  devices.  In  ongoing  work,  we  are  extending  the  model 
to  multiple  dimensions,  which  should  allow  a  more  realistic  representation  of  the  electric  field  and  better 
quantitative  comparison  with  experiment. 


Acknowledgments 

This  project  was  sponsored  in  part  by  the  Air  Force  Office  of  Scientific  Research  (monitored  by  F. 
Fahroo,  AFOSR/RSL),  and  by  a  grant  of  High  Performance  Computing  time  from  the  Air  Force  Research 
Laboratory  Major  Shared  Resource  Center.  Work  at  The  Ohio  State  University  was  funded  in  part  by 
the  Chief  Scientist  Innovative  Research  Fund  (CSIRF)  of  the  Air  Force  Research  Laboratory  Air  Vehicles 
Directorate  (AFRL/RB).  The  authors  would  like  to  thank  M.  N.  Shneider  for  helpful  discussions.  Cleared 
for  public  release,  distribution  unlimited  (88ABW-2012-6103). 


References 

1Lieberman,  M.  A.  and  Lichtenberg,  A.  J.,  Principles  of  Plasma  Discharges  and  Materials  Processing ,  J.  Wiley,  New 
York,  1994. 

2Vasilyak,  L.  M.,  Kostyuchenko,  S.  V.,  Kudryavtsev,  N.  N.,  and  Filyugin,  I.  V.,  “Fast  Ionisation  Waves  under  Electrical 
Breakdown  Conditions,”  Physics- Uspekhi,  Vol.  37,  No.  4,  1994,  pp.  247-269. 

3 Fridman,  A.  and  Kennedy,  L.  A.,  Plasma  Physics  and  Engineering ,  Taylor  and  Francis,  New  York,  2004. 

4Macheret,  S.  O.,  Shneider,  M.  N.,  and  Miles,  R.  B.,  “Modeling  of  Air  Plasma  Generation  by  Repetitive  High-Voltage 
Nanosecond  Pulses,”  IEEE  Transactions  on  Plasma  Science ,  Vol.  30,  No.  3,  2002,  pp.  1301-1314. 

5Opaits,  D.  F.,  Likhanskii,  A.  V.,  Neretti,  G.,  Zaidi,  S.,  Shneider,  M.  N.,  Miles,  R.  B.,  and  Macheret,  S.  O.,  “Experimental 
Investigation  of  Dielectric  Barrier  Discharge  Plasma  Actuators  Driven  by  Repetitive  High-Voltage  Nanosecond  Pulses  with  DC 
or  Low  Frequency  Sinusoidal  Bias,”  Journal  of  Applied  Physics,  Vol.  104,  2008,  pp.  043304. 

6Roupassov,  D.  V.,  Nikipelov,  A.  A.,  Nudnova,  M.  M.,  and  Starikovskii,  A.  Y.,  “Flow  Separation  Control  by  Plasma 
Actuator  with  Nanosecond  Pulsed-Periodic  Discharge,”  AIAA  Journal,  Vol.  47,  No.  1,  2009,  pp.  168-185. 

7(Udagawa),  K.  T.,  Zuzeek,  Y.,  Lempert,  W.  R.,  and  Adamovich,  I.  V.,  “Characterization  of  a  Surface  Dielectric  Bar¬ 
rier  Discharge  Plasma  Sustained  by  Repetitive  Nanosecond  Pulses,”  Plasma  Sources  Science  and  Technology,  Vol.  20,  2011, 
pp.  055009. 

8Little,  J.,  Takashima,  K.,  Nishihara,  M.,  Adamovich,  I.,  and  Samimy,  M.,  “Separation  Control  with  Nanosecond-Pulse- 
Driven  Dielectric  Barrier  Discharge  Plasma  Actuators,”  AIAA  Journal,  Vol.  50,  No.  2,  2012,  pp.  350-365. 

9Murray,  R.  C.,  Zaidi,  S.  H.,  Carraro,  M.  R.,  Vasilyak,  L.  M.,  Macheret,  S.  O.,  Shneider,  M.  N.,  and  Miles,  R.  B., 
“Magnetohydrodynamic  Power  Generation  Using  Externally  Ionized,  Cold,  Supersonic  Air  as  Working  Fluid,”  AIAA  Journal, 
Vol.  44,  No.  1,  2006,  pp.  119-127. 

10Adamovich,  I.  V.,  Lempert,  W.  R.,  Nishihara,  M.,  Rich,  J.  W.,  and  Utkin,  Y.  G.,  “Repetitively  Pulsed  Nonequilibrium 
Plasmas  for  Magnetohydrodynamic  Flow  Control  and  Plasma- Assisted  Combustion,”  Journal  of  Propulsion  and  Power,  Vol.  24, 
No.  6,  2008,  pp.  1198-1215. 

11Pancheshnyi,  S.  V.,  Lacoste,  D.  A.,  Bourdon,  A.,  and  Laux,  C.  O.,  “Ignition  of  Propane-Air  Mixtures  by  a  Repetitively 
Pulsed  Nanosecond  Discharge,”  IEEE  Transactions  on  Plasma  Science,  Vol.  34,  No.  6,  2006,  pp.  2478-2487. 

12Sun,  W.,  Uddi,  M.,  Won,  S.  H.,  Ombrello,  T.,  Carter,  C.,  and  Ju,  Y.,  “Kinetic  Effects  of  Non-Equilibrium  Plasma-Assisted 
Methane  Oxidation  on  Diffusion  Flame  Extinction  Limits,”  Combustion  and  Flame,  Vol.  159,  2012,  pp.  221-229. 

13Pugh,  E.  R.,  Wallace,  J.,  Jacob,  J.  H.,  Northam,  D.  B.,  and  Daugherty,  J.  D.,  “Optical  Quality  of  Pulsed  Electron-Beam 
Sustained  Lasers,”  Applied  Optics,  Vol.  13,  No.  11,  1974,  pp.  2512-2517. 

14 Aleksandrov,  V.  V.,  Koterov,  V.  N.,  Pustovalov,  V.  V.,  Soroka,  A.  M.,  and  Suchkov,  A.  F.,  “Space-Time  Evolution  of  the 
Cathode  Layer  in  Electron-Beam-Controlled  Lasers,”  Soviet  Journal  of  Quantum  Electronics,  Vol.  8,  No.  1,  1978,  pp.  59-63. 

15Popov,  N.  A.,  “Investigation  of  the  Mechanism  for  Rapid  Heating  of  Nitrogen  and  Air  in  Gas  Discharges,”  Plasma  Physics 
Reports,  Vol.  27,  No.  10,  2001,  pp.  886-896. 

16Flitti,  A.  and  Pancheshnyi,  S.,  “Gas  Heating  in  Fast  Pulsed  Discharges  in  N2-O2  Mixtures,”  The  European  Physical 
Journal:  Applied  Physics,  Vol.  45,  2009,  pp.  21001. 

17 Aleksandrov,  N.  L.,  Kindysheva,  S.  V.,  Nudnova,  M.  N.,  and  Starikovskiy,  A.  Y.,  “Mechansim  of  Ultra-Fast  Heating 
in  a  Non-Equilibrium  Weakly  Ionized  Air  Discharge  Plasma  in  High  Electric  Fields,”  Journal  of  Physics  D:  Applied  Physics, 
Vol.  43,  2010,  pp.  255201. 

18Popov,  N.,  “Fast  Gas  Heating  in  a  Nitrogen- Oxygen  Discharge  Plasma:  I.  Kinetic  Mechanism,”  Journal  of  Physics  D: 
Applied  Physics,  Vol.  44,  2011,  pp.  285201. 


9  of  16 


American  Institute  of  Aeronautics  and  Astronautics 


Downloaded  by  Igor  Adamovich  on  January  9,  2013  I  http://arc.aiaa.org  I  DOI:  10.2514/6.2013-458 


19Unfer,  T.  and  Boeuf,  J.  P.,  “Modelling  of  a  Nanosecond  Surface  Discharge  Actuator,”  Journal  of  Physics  D:  Applied 
Physics,  Vol.  42,  2009,  pp.  194017. 

20Nishihara,  M.,  Takashima,  K.,  Rich,  J.  W.,  and  Adamovich,  I.  V.,  “Mach  5  Bow  Shock  Control  by  a  Nanosecond  Pulse 
Surface  Dielectric  Barrier  Discharge,”  Physics  of  Fluids,  Vol.  23,  2011,  pp.  066101. 

21Kirpichnikov,  A.  A.  and  Starikovskii,  A.  Y.,  “Nanosecond  Pulse  Discharge  -  Always  Uniform?”  IEEE  Transactions  on 
Plasma  Science,  Vol.  36,  No.  4,  2008,  pp.  898-899. 

22Poggie,  J.,  Bisek,  N.  J.,  Adamovich,  I.  V.,  and  Nishihara,  M.,  “Numerical  Simulation  of  Nanosecond-Pulse  Electrical 
Discharges,”  AIAA  Paper  2012-1025,  American  Institute  of  Aeronautics  and  Astronautics,  Reston,  VA,  January  2012. 

23Likhanskii,  A.  V.,  Shneider,  M.  N.,  Macheret,  S.  O.,  and  Miles,  R.  B.,  “Modeling  of  Dielectric  Barrier  Discharge  Plasma 
Actuators  Driven  by  Nanosecond  Pulses,”  Physics  of  Plasmas,  Vol.  14,  2007,  pp.  073501. 

24Breden,  D.  and  Raja,  L.  L.,  “Gas  Heating  Effects  in  a  Nanosecond- Pulse  Streamer  Discharge  Interacting  with  a  Supersonic 
O2-H2  Flow,”  IEEE  Transactions  on  Plasma  Science ,  Vol.  39,  No.  11,  2011,  pp.  2250-2251. 

25 Che,  X.,  Shao,  T.,  Nie,  W.,  and  Yan,  P.,  “Numerical  Simulation  on  a  Nanosecond- Pulse  Surface  Dielectric  Barrier 
Discharge  Actuator  in  Near  Space,”  Journal  of  Physics  D:  Applied  Physics,  Vol.  45,  2012,  pp.  145201. 

26Wang,  C.-C.,  “Energy  and  Force  Prediction  for  a  Nanosecond  Pulsed  Dielectric  Barrier  Discharge  Actuator,”  Journal  of 
Applied  Physics,  Vol.  Ill,  No.  10,  2012,  pp.  103302. 

27Uddi,  M.,  Jiang,  N.,  Adamovich,  I.  V.,  and  Lempert,  W.  R.,  “Nitric  Oxide  Density  Measurements  in  Air  and  Air/Fuel 
Nanosecond  Pulse  Discharges  by  Laser  Induced  Fluorescence,”  Journal  of  Physics  D:  Applied  Physics,  Vol.  42,  2009,  pp.  075205. 

28Kossyi,  I.  A.,  Kostinsky,  A.  Y.,  Matveyev,  A.  A.,  and  Silakov,  V.  P.,  “Kinetic  Scheme  of  the  Nonequilibrium  Discharge 
in  Nitrogen- Oxygen  Mixtures,”  Plasma  Sources  Science  and  Technology ,  Vol.  1,  1992,  pp.  207-220. 

29McNeal,  R.  J.,  Jr.,  M.  E.  W.,  and  Cook,  G.  R.,  “Quenching  of  Vibrationally  Excited  N2  by  Atomic  Oxygen,”  Chemical 
Physics  Letters,  Vol.  16,  No.  3,  1972,  pp.  507-510. 

30Huxley,  L.  G.  H.  and  Crompton,  R.  W.,  The  Diffusion  and  Drift  of  Electrons  in  Gases,  Wiley,  New  York,  1974. 

31Itikawa,  Y.,  Hayashi,  M.,  Ichimura,  A.,  Onda,  K.,  Sakimoto,  K.,  Takayanagi,  K.,  Nakamura,  M.,  Nishimura,  M.,  and 
Takayanagi,  T.,  “Cross  Sections  for  Collisions  of  Electrons  and  Photons  with  Nitrogen  Molecules,”  Journal  of  Physical  and 
Chemical  Reference  Data,  Vol.  16,  1986,  pp.  985-1010. 

32Itikawa,  Y.,  Ichimura,  A.,  Onda,  K.,  Sakimoto,  K.,  Takayanagi,  K.,  Hatano,  Y.,  Hayashi,  M.,  Nishimura,  H.,  and 
Tsurubichi,  S.,  “Cross  Sections  for  Collisions  of  Electrons  and  Photons  with  Oxygen  Molecules,”  Journal  of  Physical  and 
Chemical  Reference  Data,  Vol.  18,  1989,  pp.  23-42. 

33Haydon,  S.  C.  and  Williams,  O.  M.,  “Combined  Spatial  and  Temporal  Studies  of  Ionization  Growth  in  Nitrogen,”  Journal 
of  Physics  D:  Applied  Physics,  Vol.  9,  1976,  pp.  523-536. 

34Chanrion,  O.  and  Neubert,  T.,  “A  PIC-MCC  Code  for  Simulation  of  Streamer  Propagation  in  Air,”  Journal  of  Compu¬ 
tational  Physics,  Vol.  227,  2008,  pp.  7222-7245. 

35Meyyappan,  M.  and  Kreskovsky,  J.  P.,  “Glow  Discharge  Simulation  Through  Solutions  to  the  Moments  of  the  Boltzmann 
Transport  Equation,”  Journal  of  Applied  Physics,  Vol.  68,  No.  4,  1990,  pp.  1506-1512. 

36Wilcoxson,  M.  H.  and  Manousiouthakis,  V.  I.,  “Simulation  of  a  Three-Moment  Fluid  Model  of  a  Two-Dimensional  Radio 
Frequency  Discharge,”  Chemical  Engineering  Science,  Vol.  51,  No.  7,  1996,  pp.  1089-1106. 

37Chen,  G.  and  Raja,  L.  L.,  “Fluid  Modeling  of  Electron  Heating  in  Low-Pressure,  High-Frequency  Capacitively  Coupled 
Plasma  Discharges,”  Journal  of  Applied  Physics,  Vol.  96,  No.  11,  2004,  pp.  6073-6081. 

38 Adamovich,  I.  V.  and  Rich,  J.  W.,  “Emission  and  Shock  Visualization  in  Nonequilibrium  Nitrogen  Afterglow  Plasma,” 
Journal  of  Applied  Physics,  Vol.  102,  2007,  pp.  083303. 

39Hakim,  A.  and  Shumlak,  U.,  “Two-Fluid  Physics  and  Field- Reversed  Configurations,”  Physics  of  Plasmas,  Vol.  14,  2007, 
pp.  055911. 

40Poggie,  J.,  Adamovich,  I.,  Bisek,  N.,  and  Nishihara,  M.,  “Numerical  Simulation  of  Nanosecond-Pulse  Electrical  Dis¬ 
charges,”  Plasma  Sources  Science  and  Technology ,  Vol.  22,  2012,  pp.  015001. 

41Burgers,  J.  M.,  Flow  Equations  for  Composite  Gases,  Academic  Press,  New  York,  1969. 

42Golant,  V.  E.,  Zhilinsky,  A.  P.,  and  Sakharov,  I.  E.,  Fundamentals  of  Plasma  Physics,  J.  Wiley,  New  York,  1980,  Trans. 
K.  Z.  Vedeneyeva  and  V.  F.  Agranat. 

43 Williams,  F.  A.,  Combustion  Theory,  Addison- Wesley,  Menlo  Park  CA,  4th  ed.,  1988. 

44White,  F.  M.,  Viscous  Fluid  Flow,  McGraw-Hill,  New  York,  2nd  ed.,  1991. 

45Mahadevan,  S.  and  Raja,  L.,  “Simulations  of  Direct- Current  Air  Glow  Discharges  at  Pressures  ~  1  Torr:  Discharge 
Model  Validation,”  Journal  of  Applied  Physics,  Vol.  107,  2010,  pp.  093304. 

46 Chen,  N.  H.  and  Othmer,  D.  F.,  “New  Generalized  Equation  for  Gas  Diffusion  Coefficient,”  Journal  of  Chemical  and 
Engineering  Data,  Vol.  7,  No.  1,  1962,  pp.  37-41. 

47Viehland,  L.  A.  and  Mason,  E.  A.,  “Transport  Properties  of  Gaseous  Ions  over  a  Wide  Energy  Range,  IV,”  Atomic  Data 
and  Nuclear  Data  Tables,  Vol.  60,  1995,  pp.  37-95. 

48Nelson,  D.,  Benhenni,  M.,  Eichwald,  O.,  and  Yousfi,  M.,  “Ion  Swarm  Data  for  Electrical  Discharge  Modeling  in  Air  and 
Flue  Gas  Mixtures,”  Journal  of  Applied  Physics,  Vol.  94,  No.  1,  2003,  pp.  96-103. 

49Deconinck,  T.,  Mahadevan,  S.,  and  Raja,  L.  L.,  “Simulation  of  Direct-Current  Surface  Plasma  Discharge  Phenomena  in 
High-Speed  Flow  Actuation,”  IEEE  Transactions  on  Plasma  Science ,  Vol.  35,  No.  5,  2007,  pp.  1301-1311. 

50Ivanov,  A.  V.,  Takhtenberg,  S.,  Bertram,  A.  K.,  Gershenzon,  Y.  M.,  and  Molina,  M.  J.,  “OH,  HO2,  and  Ozone  Gaseous 
Diffusion  Coefficients,”  Journal  of  Physical  Chemistry  A,  Vol.  Ill,  2007,  pp.  1632-1637. 

51Breden,  D.  and  Raja,  L.  L.,  “Simulations  of  Nanosecond  Pulsed  Plasmas  in  Supersonic  Flow,”  AIAA  Paper  2009-3594, 
American  Institute  of  Aeronautics  and  Astronautics,  Reston,  VA,  June  2009. 

52Poggie,  J.,  “Numerical  Simulation  of  Direct  Current  Glow  Discharges  for  High-Speed  Flow  Control,”  Journal  of  Propulsion 
and  Power,  Vol.  24,  No.  5,  2008,  pp.  916-922. 


10  of  16 


American  Institute  of  Aeronautics  and  Astronautics 


Downloaded  by  Igor  Adamovich  on  January  9,  2013  I  http://arc.aiaa.org  I  DOI:  10.2514/6.2013-458 


53Poggie,  J.,  “High-Order  Numerical  Methods  for  Electrical  Discharge  Modeling,”  AIAA  Paper  2010-4632,  American  In¬ 
stitute  of  Aeronautics  and  Astronautics,  Reston  VA,  2011. 

54Beam,  R.  and  Warming,  R.,  “An  Implicit  Factored  Scheme  for  the  Compressible  Navier-Stokes  Equations,”  AIAA 
Journal ,  Vol.  16,  No.  4,  1978,  pp.  393-402. 

55 Morrison,  J.,  “Flux  Difference  Split  Scheme  for  Turbulent  Transport  Equations,”  AIAA  Paper  90-5251,  American  Institute 
of  Aeronautics  and  Astronautics,  Reston,  VA,  October  1990. 

56Walters,  R.  W.,  Cinnella,  P.,  Slack,  D.  C.,  and  Halt,  D.,  “Characteristic-Based  Algorithms  for  Flows  in  Thermo-Chemical 
Nonequilibrium,”  AIAA  Paper  90-0393,  American  Institute  of  Aeronautics  and  Astronautics,  Reston,  VA,  January  1990. 

57Gaitonde,  D.  and  Shang,  J.  S.,  “Accuracy  of  Flux-Split  Algorithms  in  High-Speed  Viscous  Flows,”  AIAA  Journal ,  Vol.  31, 
No.  7,  1993,  pp.  1215-1221. 

58  Anderson,  W.  K.,  Thomas,  J.  L.,  and  van  Leer,  B.,  “A  Comparison  of  Finite  Volume  Flux  Vector  Splittings  for  the  Euler 
Equations,”  AIAA  Paper  85-0122,  American  Institute  of  Aeronautics  and  Astronautics,  Reston,  VA,  January  1985. 

59Ventzek,  P.  L.  G.,  Hoekstra,  R.  J.,  and  Kushner,  M.  J.,  “Two-Dimensional  Modeling  of  High  Plasma  Density  Inductively 
Coupled  Sources  for  Materials  Processing,”  Journal  of  Vacuum  Science  and  Technology  B ,  Vol.  12,  No.  1,  1994,  pp.  461-477. 

60Holst,  T.  L.,  “Transonic  Flow  Computations  Using  Nonlinear  Potential  Methods,”  Progress  in  Aerospace  Sciences , 
Vol.  36,  2000,  pp.  1-61. 

61  Soloviev,  V.  R.  and  Krivtsov,  V.  M.,  “Surface  Barrier  Discharge  Modelling  for  Aerodynamic  Applications,”  Journal  of 
Physics  D:  Applied  Physics ,  Vol.  42,  2009,  pp.  125208. 


No. 

Species 

'I/O 

rLs 

(eV /particle) 

0’s 

(10-20  m2) 

1 

n2 

0 

41 

2 

o2 

0 

39 

3 

o 

2.58 

32 

4 

o3 

1.48 

45 

5 

NO 

0.94 

31 

6 

N 

4.90 

38 

7 

OfD) 

4.55 

32 

8 

N2(A3E) 

6.17 

41 

9 

N2(B3n) 

7.35 

41 

10 

N2(a,1E) 

8.40 

41 

11 

N2(C3n) 

11.03 

41 

12 

e_ 

0 

- 

13 

Nj 

16.37 

148 

14 

o+ 

16.26 

80 

15 

02+ 

12.69 

115 

16 

N2(u  =  1) 

0.29 

41 

17 

N2(v  =  2) 

0.57 

41 

18 

N2(v  =  3) 

0.86 

41 

19 

N2(-u  =  4) 

1.14 

41 

20 

N2(u  =  5) 

1.41 

41 

21 

N2(d  =  6) 

1.68 

41 

22 

N2(u  =  7) 

1.95 

41 

23 

N2(-u  =  8) 

2.21 

41 

Table  1.  Species  included  in  kinetic  model.  Heat  of  formation  referenced  to  300  K. 
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Table  2.  Reaction  mechanism.  Units  consistent  with  number  densities  in  cm  3,  and  temperatures  in  K, 
one-body  rates  in  s— 1,  two-body  rates  in  cm3/s,  and  three-body  rates  in  cm6/s. 
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12.69 +  3fcsTe/(2e) 

3 

O2  +  e_  —¥  0  +  0+  +  e_  +  e_ 

18.84 +  3fcsTe/(2e) 

4 

N2  +  e"  -4  N2(A3E)  +  e" 

6.17 

5 

N2  +  e“  -4  N2(B3n)  +  e" 

7.35 

6 

N2  +  e“  — >  N2(C3n)  +  e" 

11.03 

7 

N2  +  e-  N2(a,1E)  +  e“ 

8.40 

8 

N2  +  e“  ->  N  +  N  +  e_ 

9.80 

9 

02  +  e_  — >  0  +  0  +  e“ 

5.17 

10 

02  +  e“  ->  0  +  0(1D)  +  e“ 

7.14 

43 

N2(v  =  0)  +  e_  — >  N2(v  =  1)  +  e_ 

0.29 

44 

N2(v  =  0)  +  e_  — >  N2(v  =  2)  +  e_ 

0.57 

45 

N2(v  =  0)  +  e_  — >  N2(v  =  3)  +  e_ 

0.86 

46 

N2(v  =  0)  +  e_  — >  N2(v  =  4)  +  e_ 

1.14 

47 

N2(v  =  0)  +  e_  — >  N2(v  =  5)  +  e_ 

1.41 

48 

N2(v  =  0)  +  e_  — >  N2(v  =  6)  +  e_ 

1.68 

49 

N2(v  =  0)  +  e_  — >  N2(v  =  7)  +  e~ 

1.95 

50 

N2(v  =  0)  +  e  — >  N2(v  =  8)  +  e 

2.21 

Table  3.  Energy  lost  by  electrons  in  inelastic  collisions. 
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(a)  Rates  of  ionization  and  electronic  excitation  (Reac¬ 
tions  1-10). 


(b)  Rates  of  vibrational  excitation  (Reactions  43-50). 


(c)  Electron- neutral  collision  rate. 


Figure  2.  Test  of  Boltzmann  solver:  predicted 
Townsend  ionization  coefficient  for  nitrogen  compared 
to  experimental  data  of  Haydon  and  Williams.33 


Figure  3.  Diagram  of  computational  domain. 


Figure  1.  Quantities  determined  from  zero¬ 
dimensional  solution  of  Boltzmann  equation. 
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(c)  Temperatures  (/c^Ts/e). 

Figure  4.  Discharge  profiles  at  12  ns  for  Vo  =  5  kV. 
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(c)  Temperatures  (/c#Ts/e). 

Figure  5.  Discharge  profiles  at  24  ns  for  Vo  = 
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Electron  Velocity  (m/s) 
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(a)  Temperature,  Vq  =  5  kV  case. 


(b)  Velocity,  Vq  =  5  kV  case. 
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(c)  Temperature,  Vo  =  27  kV  case. 


(d)  Velocity,  Vo  =  27  kV  case. 


Figure  6.  Profiles  of  the  properties  of  the  bulk  gas  for  selected  times  in  the  simulation.  Arrows  indicate 
direction  of  wave  motion. 
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